Quantum effective potential, electron transport and conformons in biopolymers 
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In the Kirchhoff model of a biopolymer, conformation dynamics can be described in terms of 
solitary waves, for certain special cross-section asymmetries. Applying this to the problem of electron 
transport, we show that the quantum effective potential arising due to the bends and twists of the 
polymer enables us to formalize and quantify the concept of a conformon that has been hypothesized 
in biology. Its connection to the soliton solution of the cubic nonlinear Schrodinger equation emerges 
in a natural fashion. 
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Geometry and topology of long chain biopolymers such as proteins and DNA play a significant role |l| during 
processes such as replication and transcription. Interesting experiments for studying the conformation and elastic 
properties of a single polymer by bending or twisting it have been devised [2J . The static properties of semi-flexible 
biopolymers such as actin which have only bending energy, and are described by the well-known wormlike chain 
0^ ' (WLC) model with a single elastic constant, the bending modulus. In contrast, the static DNA with its double- 
helix structure is described by the wormlike rod chain (WLRC) model Q| with an additional elastic constant, the twist 
rigidity. Although considerable work has been done on various equilibrium properties of both these elastic models, 
their intrinsic dynamical properties have not been studied so far. The latter play a crucial role in the mechanisms of 
' energy and information propagation along a biopolymer ,an issue of vital interest to biologists, chemists and physicists 
r**"" , alike. The equilibrium properties have been studied in U The study of intrinsic static and dynamical properties 
of biopolymers, taking into account their geometry, is a subject of great importance Such issues are of vital 

interest in biology as well as physics, since they would help us understand the mechanisms of storage and transport 
of energy and charge along a biopolymer. 

In this paper, we describe a biopolymer using the Kirchhoff model 5]. This model starts with equations that govern 
the dynamics of a thin rod that in fact characterizes a polymer in the well established WLRC model mentioned above. 
We are interested in its intrinsic dynamics, as well as its effect on electron transport, since the measured electrical 
conductivities of certain polymers are seen to be much larger than expected due to conventional mechanisms. 7] Under 
certain conditions, polymer conformations take on the form of spatially localized nonlinear excitations. Applying this 
to the problem of electron transport, we show that the quantum effective potential arising due to the bends and twists 
of the polymer enables us to formalize and quantify the concept of a conformon that has been put forward in biology 
0- 0- It is expected to play an important role in statics and dynamics of biopolymers in general. Its connection to 
. i-H , the soliton solution of the cubic nonlinear Schrodinger equation emerges in a natural fashion. 

We consider the biopolymer to be a very thin elastic filament (or rod) modeled by a sinpjg, which is defined as 
a space curve R(s,i), along with a smooth unit vector field d 2 , perpendicular to the curve. Here s denotes the arc 
length of the polymer and t is the time. The unit tangent to the curve is given by d 3 and the third unit vector of 
the triad is defined as di = d 2 x d 3 , so that the triad (d.3,d.2,di) forms a right-handed, orthonormal frame at every 
point on the curve. 

The space derivatives of the vectors of the frame can be shown to be given by the compact expression 

di, s =kxd„ (1) 

where i = 1, 2, 3, the subscript s stands for 4-, and k, the Darboux vector or "twist" vector is given by 

k(s,t) = fcidi +k 2 d 2 + k 3 d 3 . (2) 

Its components fc,-, i — 1, 2, 3, can be expressed as a function of the curvature k, torsion r and the angle <f> between 
the principal normal to the curve and di, the normal to the strip: 

(fei, ki^ks) = (ksha.(f), kcoscf), t + <j> 3 ) (3) 

Here, 



fc = |d 3 , s | 



(4) 
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and 

r = d 3 .(d 3:S x d 3 , ss )/fc 2 . (5) 

The Kirchhoff equations that govern the dynamics of the biopolymer (modelled as a thin elastic rod) are given (in 
their dimensionless form) by 

g s = Rtt (6) 

and 

m s + d 3 x g = adi x d Mt + d 2 x d 2 ,tt, (7) 

with 

m = fcidi + ak 2 d 2 + bk 3 d 3 , (8) 

where the subscript t stands for the time derivative ^. In these equations g(s, t) and m(s, t) represent the force (or 
tension) and the torque acting on each cross-section of the rod. The equations are obtained from the conservation 
of linear and angular momentum. The parameter a(0<a<l)isa measure of the bending asymmetry of its 
cross-section, b = 2a/(l + cr)(l + a), a being the Poisson ratio, is a measure of the change in volume of the rod as it 
is stretched. 

First we consider the static version of the Kirchhoff equations Using the general expression 

g = gidi + g 2 d 2 + g 3 d 3 (9) 
leads to the following system of equations [l(i|: 

9i,s + k 2 g 3 - k 3 g 2 = (10) 

.92, s + hgi - kig 3 = (11) 

53, s + kig 2 - k 2 gi = (12) 

.92 = fci, s + {b- a)k 2 k 3 (13) 

gi = -ak 2 , s + (b - l)feife 3 (14) 

6fc 3 , s + (a- l)hk 2 = (15) 

For all 4> — mr/2, n an integer, using eq. in eq. 115|l shows that 

k 3 = t = t q . (16) 

Thus the torsion of the polymer is a constant, denoted by tq. As an example, we first take <j) — -k in eq. @ and 
analyze eqs. (I1UI15|) . We find b — 2a. This implies a = —1/(1 + er). Further, 

g = a fc s di + (a - b) t kd 2 + a(-^k 2 + C 2 )d 3 , (17) 

where C 2 is an integration constant. To understand its physical significance, note that for k = (a straight polymer), 
g = a C 2 d 3 . This essentially means that C 2 represents the tension in the polymer. 
With this result, eq. IjlOfl leads to the following equation for the curvature k: 
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k ss + y = (C 2 - r 2 )fc (18) 

Equation l|18|ll has two trivial solutions: the straight line k = 0, and the circular helix k = ^/2{C 2 — Tq)- More 
interestingly, it admits the following nontrivial solution: 



k = 2^Jc 2 -t 2 sech yfc 2 -T* s, (19) 

where as already stated, To and C 2 are constants. 

For 4> = 7^ , following the same procedure, we get 6 = 2. This implies a = —(1 + cr)/a. We can show that k satisfies 
an equation of the same form as eq. I|18fl . In fact, we can verify that for all <j> = mr/2, n any integer, the curvature has 
the form given in eq. I|19|l . and as already found, r is just a constant, r . Since the curvature must be real, eq. I|19|) 
shows that C2 must be always greater than t 2 . In the case of a planar polymer, To = 0. Thus physically, a larger 
tension is needed to get a twisted polymer, for the case under discussion. 

Turning our attention to dynamical solutions, [f| have noted that the Kirchhoff equations 10 and Q can support 
traveling wave solutions for the curvature k, called Kovalevskaya waves. These are of the same form as the static 
solution 1)19(1 . where now s is replaced by £ = (s — vt), with v the speed of these spatially localized, solitary waves, 
which propagate without change of form. These arise due to a certain nontrivial scaling property || satisfied by 
Kirchhoff equations. 

We will consider possible quantum mechanical implications of this non-trivial solution for k, with regard to electron 
transport on a biopolymer. It has been shown by jllj . [13 and that a quantum particle in a thin tube whose axis 
follows a space curve with curvature k and constant torsion To (as in our case) feels an effective potential 0] of the 
form 



"<"<■>- 51-^41- (20 > 

Writing down the Schrodinger equation for an electron in the presence of the above effective potential, and mak- 
ing a gauge transformation of the wave function ipi, by using the following appropriate phase factor ipi(s,t) = 

2 

^(s, t) exp(— ihj^t), we obtain 

h 2 ( d 2 k 2 (s)\ ( . . d , 

' » ip(s,t)=in—ip(s,t) (21) 



2m \ds 2 4 / rv ' ' dt 
After rescaling the time such that -£-t — > u and the coordinate s — * y/2si, the Schrodinger equation reads: 

k 2 

i 4>« + *Ps lsl + = 0, (22) 

where k — fe(si), and the subscripts si and u stand for the partial derivatives and J^. 
Looking for solutions of eq. (|22|) of the form 

ip(si, u) = k(si) expect u), (23) 

we get 

k SlSl +^pj=ak. (24) 

This equation has the same form as eq. (|18|) . provided 

a=(C 2 -r 2 ). (25) 
Using the solution given in eq. (|19|l (with s replaced bysi) in eq. ll'.'il) . we get 

ip(s±,u) = 2ya sech \fa s\ exp i(a u), (26) 



FIG. 1: Polymer conformation for curvature k as given in eq. 119H with (C2 — tq ) = a = 1 and torsion r = To = 1. Note the 
localized twisted loop on the polymer. 



where a > 0. It is readily seen that the wave function of the electron is localized around that point on the polymer 
where the maximum of its curvature is located. Further, it has a simple sinusoidal time-dependence like a " breathed . 
Since a = k 2 (s = 0)/4 = fc§/4, eq. leads to 

f+-o 2 = C 2 . (27) 

This leads to an interesting constraint between the maximum curvature ko and the constant torsion to of the polymer, 
C2 being the constant representing tension. 

We choose various values of a and To that satisfy this constraint. Note that a < C2, from eq. I|25|l. The actual 
conformation of the polymer which has a space-dependent curvature k — 2 v / a sech y/a Si and a constant torsion 
t = To, can be found by integrating eqs. ijljl. Typically, we find that polymer has a single non-intersecting twisted 
loop, centered around si = 0. It straightens out as si — > ±00, as it should, since its curvature is readily seen to 
vanish in those limits. Figure 1 gives an example of such a conformation, for Ci = 2, with a = 1 and To = 1. For the 
same C2, smaller values of torsion, e.g., To = 0.7, make the loop curve more around the center, while for larger values, 
e.g., t = 1.23, the opposite happens, and the loop starts "unraveling" and straightens out more. Our results show 
how the above conformation of a polymer that emerges directly from static Kirchhoff equations, can lead to electron 
localization, i.e., "trapping" of an electron around the maximum curvature point on the twisted loop that develops 
mid-way on the polymer. As already mentioned, the dynamical solutions for the curvature k are just Kovalevskaya 
traveling waves, given by 

k(si,u) = k(si — v u) = 2\/a sech[vcv (si — v u)]. (28) 

For this case, the wave function of the electron is to be found as the solution of the corresponding time-dependent 
Schrodinger equation eq. (|22[). where now k — fc(sj — v u) is given in l|28|) . Thus for this dynamical case, we look for 
a solution of the form 

ip(si,u) = k(si — v u) exp i[X si + fx u]. (29) 
where A and /1 are to be found by substituting eq. (|29() into eq. I|22|l . with k as in l|28() . After some algebra, we find 

\ = (v/2) ; M =( a -A 2 ) = [C 2 -T 2 -^], (30) 
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on using eq. Q25[l. Substituting for A and /j, from eq. I|30[) and k(s% — vu) from eq. Q28[l . eq. (|29|l becomes 

2 

V-> = sech[-\fa{s\ — vu)] expi[— si + (a — (31) 

This is identical to the envelope soliton solution of the following, completely integrable cubic nonlinear Schrodinger 
equation (CNLSE), 

\ib\ 2 

iipu+ip SlSl + ^^ = 0- (32) 

This is as expected, because for the solution (|29H that we have considered, k 2 = \ip\ 2 , so that eq. (|22|l reduces to 
eq. 

^From eq. (|31|) . it is clear that the envelope soliton has a localized profile: Its modulus travels with envelope velocity 
V e = v, while its phase has a carrier velocity V c given by 

v.-M-'ZiM- r-**-™ ,33, 

2v 2v 

This leads to the well known inequality v[v — 2V C ] > between these two velocities of the CNLSE soliton on the 
polymer. From eq. (|33(l . we see that for a given v, V c depends on the tension and torsion of the polymer. 

Incorporating the additional phase factor exp(— it 2 u) due to the gauge transformation we had made earlier, we 
finally obtain the following travelling wave solution for the wave function ipi(s\,u): 

ipi(si,u) — lyfa. sech {\fa (si — vu)) 

exp i[v si + (C 2 - 2t 2 - v 2 )u], (34) 

where a is defined in eq. I|25|l . It is easy to see that in this case, the electron gets trapped by a moving potential 
well, which travels along the polymer. To understand the conformation here, we note that the polymer now has a 
curvature which is a Kovalevskaya solitary wave, traveling without change of form: k = 2y / asech(y / a< (si — vu). The 
conformation is again a twisted loop, but now it travels with a constant velocity v. Thus the electron gets trapped in 
the loop, and is transported along with it, on the polymer. As we have seen, its transport is soliton-like in this case. 

We believe that our results provide a precise dynamical underpinning for the conformon concept hypothesized by 
various authors 0, to play an important role in biology. Green and Ji § state that a conformon is a localized 

packet of energy (and genetic information). It is an energy packet associated with a conformational strain, which 
is localized in a region much shorter than the length of the molecule Q. We find the curvature A: to be a localized 
function. Since the energy density on the polymer is proportional to fc 2 , this leads to a localized packet of energy. 

Volkenstein|l5j suggests that a conformon is like an "electron plus conformational change". Kemeny and Goklany 
remark that "in some sense, the conformon is a generalization of a polaron". As is well known, a polaron is a 
localized electronic bound state in a discrete lattice, which is not perfectly periodic. It is formed by the trapping of 
the electron due to the nonlinearities arising from its strong coupling to the lattice (phonons). Here, we mention 
that in the specific context of an a-helical protein, starting with a quantum mechanical discrete lattice model, and 
invoking electron-phonon coupling, a CNLSE has been derived in the continuum approximation, by Davvdov|17|. 

On the other hand, our work deals with conformational aspects. Using the Kirchhoff model (which is a continuum 
model per se), we have shown that a localized electronic state arises in the curved and twisted polymer. This is 
essentially because its curvature and torsion "interact" with the electron by inducing a potential well, which traps it, 
in addition to creating a nonlinearity in its Schrodinger equation. While this scenario is indeed somewhat analogous 
to the polaron picture described above, the origin of the two mechanisms are quite distinct, with the curved geometry 
of the polymer playing a key role in the creation of a conformon. We conjecture that the moving soliton solution that 
arises, along with its robust propagation can provide an explanation for the unexpectedly high electrical conductivity 
(around 10 22 mho/cm) found in certain biopolymers. It represents a novel mechanism of charge transport without 
dissipation which is not restricted to low temperatures. 

Finally, we have shown how the geometry of polymers, nonlinearity and quantum particle transport are intimately 
related. Thus our results are also likely to be of significance in other kinds of transport phenomena in molecular biology. 
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